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Quantum states in the Earth's gravitational field were observed, when ultra-cold neutrons fall un- 
' der gravity. The experimental results can be described by the quantum mechanical scattering model 

as it is presented here. We also discuss other geometries of the experimental setup which correspond 
to the absence or the reversion of gravity. Since our quantum mechanical model describes, particu- 
larly, the experimentally realized situation of reversed gravity quantitatively, we can practically rule 
out alternative explanations of the quantum states in terms of pure confinement effects. 

Q • PACS numbers: 03.65. Ge,03.65.Ta,04.62.+v,04.80.-y,61.12.Ex 

Si , 

I. INTRODUCTION 

X 

A gravitationally bound quantum system has been realized experimentally with ultra-cold neutrons 
falling under gravity and reflecting off a "neutron mirror" QQ. UCN are neutrons which, in contrast to 
faster neutrons, are reflected at all angles of incidence. For such UCN, flat surfaces thus act as mirrors. 
Using an efficient neutron absorber for the removal of higher unwanted states, only neutrons in selected 
energy states are taken. This idea of observing quantum effects occurring when ultra-cold neutrons 
are stored on a plane matter surface was first discussed by V.I. Lushikov and LA. Frank 3] with the 
first concrete experimental realization in Q . An experiment in some aspects similar was discussed by H. 
Wallis et al. p| in the context of trapping atoms in a gravitational cavity. The toy model of a Schrodinger 
quantum particle bouncing in a linear gravitational field is known as the quantum bouncer 0, ]t| • 
Retroreflectors for atoms have used the electric dipole force in an evanescent light wave 0, or they 
are based on the gradient of the magnetic dipole interaction, which has the advantage of not requiring 
a laser [loj . 

A unique side-effect of the experiment with neutrons is its sensitivity to gravity-like forces at length 
scales below 10 fxm while all electromagnetic effects are extremely suppressed [ill H^. The quantum 
states probe Newtonian gravity between 10 -9 and 10 -5 m and the experiment places limits for gravity- 
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like forces there. In light of recent theoretical developments in higher-dimensional field theory (see 
also 14\ for explicit realizations in string theory), gauge fields could mediate forces that are 10 6 to 
10 12 times stronger than gravity at submillimeter distances, exactly in the interesting range of this 
experiment and might give a signal in an improved setup. 

In this article, we provide the details of a quantum mechanical calculation |15| for our experiment, 
where gravitational bound quantum states are observed for the first time. The experiment consists 
of a reflector for neutrons, called neutron mirror, an absorber for unwanted neutrons and a neutron 
detector. In our previous papers 0, 0, El the experiment and a first treatment of the data were 
presented. Fundamental limits for the spatial resolution and a first ansatz to incorporate the neutron 
scatterer can be found in . In another work |T^ a description of the neutron loss from first principles 
was developed where the rough edges of the absorber surface are treated as a time-dependent variation of 
a flat absorber position, modeling the neutron loss mechanism as a process equivalent to the ionization 
of a particle, initially confined in a well with an oscillating wall. Within the older and more simple 
model we present in this paper we are able to describe the experimental data with one micro-physical 
fit parameter which parameterizes the micro-physics of the neutron scatterer/absorber. At the moment 
our model is the only one yet in which both the data in Fig. 2 and Fig. 4 are described. 

The layout of this paper is the following: After describing the experiment and the observations, 
we recall the quantum mechanics of gravitational bound states on a free mirror, that is without an 
absorber/scatterer. Then we present an approach to describe the state selection by deriving a neutron 
loss rate due to non-specular scattering from the rough surface of our absorber. This approach explains 
the non-classical dependence of the transmission of the mirror-absorber system as a function of the 
height I of the absorber above the mirror, see Fig. 2] and Fig. [3 We take into account the deformation 
of the bound state wave functions due to the matter bulk of the state selector, and compare the full 
prediction with actual data. We also show that the data are only understood when gravity is present. 




FIG. 1: Schematic view with mirrors, absorber and 
quantum mechanical boundary conditions. In the ex- 
periment, one mirror of length 10 cm or, as an option 
as shown here, two bottom mirrors of length 6 cm were 
used. 



FIG. 2: Circles: Data from the 2 nd run 2002 with 
one bottom mirror 

El- 

Solid: Transmission coefficient 
from the phenomenological scattering model. Dash: 
The classical expectation for the neutron transmission 
coefficient. 



II. OBSERVATION OF QUANTUM STATES AND SETUP 



A description of the experiment at the Institute Laue-Langevin (ILL) can be found in It was 
installed at the UCN facility PF2 of the Institute. Here, neutrons have a velocity of several meters 
per second. They are then guided to the experiment via a curved neutron guide with a diameter of 80 
mm. At the entrance of the experiment, a collimation system cuts down to an adjustable transversal 
energy E± in the pico-eV range. Either one solid block with dimensions 10 cm x 10 cm x 3 cm or 
two solid blocks with dimensions 10 cm x 6 cm x 3 cm composed of optical glass serve as mirrors for 
UCN neutron reflection. To select different states an absorber/scatterer, a rough mirror coated with an 
Neutron absorbing alloy, is placed above the first mirror. We can vary the height I above the mirror, 
that is the size of the slit. The collimation system in front of the state selector is adjusted in such a way 
that neutrons on classical trajectories entering the experiment have to hit the mirror surface at least two 
times. After the second mirror we placed a 3 He counter for neutron detection. Fig. 1 shows a schematic 
view of our setup. Signatures of quantum states are observed in the following way: The 3 He counter 
measures the total neutron transmission F, when neutrons are traversing the mirror absorber-system 



as described. The transmission is measured as a function of the absorber height I and thus as a function 
of neutron energy since the height acts as a selector for the energy E± of the vertical motion. From 
the classical point of view, the transmission F of neutrons is proportional to the phase space volume 
allowed by the absorber. It is governed by a power law F ~ l n and n = 3/2. 

The measurements show the following: Above an Absorber/Scatterer height of about 60 /im, the 
measured transmission is in agreement with the classical expectation but below 50 /im, a deviation is 
clearly visible. Below about 15 /im, no neutrons can pass the slit. In the next section we will find that 
the vertical extension of the gravitational bound states increases with the quantum number. Ideally, we 
expect a stepwise dependence of F as a function of I. If I is smaller than the spatial width of the lowest 
quantum state, then F will be zero. When I is equal to the spatial width of the lowest quantum state 
F will increase sharply. A further increase in I should not increase F as long as I is smaller than the 
spatial width of the second quantum state. Then again, F should increase stepwise. At sufficiently high 
slit width one approaches the classical dependence. Fig. 2 shows details of the quantum regime below 
an absorber height of I — 50 /im. The transmission function depends on the horizontal neutron velocity 
and the absorption efficiency. It was found, that except for the ground state, the stepwise increase is 
mostly washed out. 



III. QUANTUM MECHANICAL DESCRIPTION OF GRAVITATION ALLY BOUND 

STATES 



The quantum mechanical treatment of a reflecting neutron mirror, made from glass, is simple. The 
glass is described by a Fermi pseudo-potential (V — iW). This potential is essentially real (\W\ << \V\) 
because of the small absorption cross section of glass and V = 100 neV is large compared with the 
transversal energy E± of the neutrons. Therefore, the potential V is set to infinity at height z = 0. 
Neutrons which hit the glass surface undergo specular reflections. 

We start with the description of the free states. On a perfect mirror, no mixing of momentum 
components take place, which leads to a decoupled one-dimensional stationary Schrodinger equation, 



2m A + v( z ))* n = En * n ,v( z) = y^ ; — (1) 

with wave functions ^ n for energies E n and the potential V(z). m is the mass of the neutron and g is 
the acceleration in the earth's gravitational field. It is convenient to use 

C = — and, above the mirror, V = mgRC,. (2) 
R 



Here, R is a scaling factor, defined as 



1/3 

R = ( ) ■ ( 3 ) 



Solutions ^ n ,g (r, t) of Eq. are obtained with an Airy function 

V n , g (r,t) =<j>{x,y) W 9 (C) e- 1 ^ 
VV S (0 = M (C - Cn) (4) 

The displacement Cn of the n-th cigcnfunction has to coincide with the n-th zero of the Airy function 
(Ai(-Cn)=0) to fulfill the boundary condition ^ n (0)=0 at the mirror. The corresponding energies E n 
with z n = i?Cn are 

E n = mgz n . (5) 
In the WKB-approximation we have to leading order 



3tt / 1 



2 / 3 

which coincides with the exact eigenvalues to better than 1 % even for the ground state 0|. The z n 
correspond to the highest point of a classical neutron trajectory with energy E n . For example, energies 
of the lowest levels (n = 1, 2, 3, 4) are 1.44 peV, 2.53 peV, 3.42 peV and 4.21 peV. The corresponding 
classical turning points 13.7 /im, 24.1 /im, 32.5 [im and 40.1 \xm. 



The aim of this experiment was to populate only some of the lowest allowed gravitationally bound 
quantum states. Higher states were removed with the absorber/scatterer at a certain height /. 



IV. PHENOMENOLOGICAL SCATTERING MODEL OF NEUTRON LOSS 



It will be most convenient to start with a polychromatic neutron beam of (locally) plane waves 
entering the system. As is well known, gaussian wave packets being closer to a particle view of a neutron 
can easily be decomposed into a Fourier integral over plane waves. Using all boundary conditions given 
by Fig. 1, one arrives at a set of matching conditions 0|. The neutron transmission of the system 
depends both on the eigenvalues and the matching conditions. Furthermore, if two bottom mirrors are 
used and shifted relative to each other a few /im in height (as it was in the 1999 beam time), there is 
an additional boundary that would change the population of the eigenstates. 

The insertion of the states yields a system of linear equations for the matching constants and its 
solution yields finally the transmission coefficient of the n th final bound state of region III. The initial 
population of the bound states of the wave guide system at coupling-out is uniform, if the vertical 
velocity distribution of the arriving beam is sufficiently wide and flat - which is of course true in our 
case, where we have about 20 cm/s or 50 peV spread in the vertical components of the arriving beam 
compared to a few peV for the lowest vertically bound states inside the mirror-absorber system [15| (see 
also [13). 

In addition, we have to introduce repopulation coefficients pj which allow us to take into account 
a step between region Ha and lib into account. If there is no step, all pj =1. In the 1999 beam time, 
the 2 nd mirror has been shifted downwards by 5 /im relative to the 1 st mirror. The matrix of overlap 
integrals of wave functions at the edge Ila/IIb is sufficiently diagonal that we can neglect the off-diagonal 
elements. We can set p\ = 0.25 and pj = 1 for j > 1. Hence, the relative shift of the bottom mirrors 
offers a possibility of controlling the relative population of in particular the ground state in the earth's 
gravitational field. However, even in the setup without step, we find a reduction of the ground state for 
unknown reasons and keep p\ as a free parameter. 

In the following, the roughness is described as an additional loss channel in the one-dimensional 
Schrodinger equation. The neutron loss is then understood in terms of non-specular scattering of 
neutrons into highly excited states which due to their large vertical energy are rapidly lost inside 
the glass of the mirror and the absorber/scatterer. Scattering happens at the rough and (due to 
Kibsorbcr — 10~ 8 eV E„ for states with low n) hard surface of the absorber, which notion shall 
enable us to derive the scattering-induced loss rate with just one undetermined micro-physical quantity 
which would be given in a micro-physical scattering model as a function of the scattering cross section. 
This parameter will be determined in the end by a fit to the data. 

The deformation of the wave functions compared with the purely gravitationally bound states due 
to the large real part of the absorber potential leads to an approximate vanishing of the bound states at 
the absorber surface. Therefore, the loss rate is calculated in terms of the deformed states. In a second 
step, one might model the neutron loss in terms of the physical parameters of the absorber surface but 
this is outside of the scope of this article. 

The neutron removal processes are modeled as a general phenomenological loss rate T n (l) of the 
n th bound state, which is taken to be proportional to the probability density of the neutrons at the 
absorber/scatterer. Here I again denotes the position of the absorber/scatterer above the mirror. The 
modulus of the bound state is then no longer constant in time since it is given as the solution at 
first order of a differential equation that determines the change of the norm to be proportional to its 
momentarily value as well as the loss rate 

d(lp n \lpn) = -(V'nIV'n) ' T n (l) ■ dt , (7) 

which yields 

GM^«> = \Pn(t)\ 2 , P n {t) = e-i r ™«-* . (8) 

The roughness which causes the loss due to scattering can be thought of as being confined to a 
region of about 2 • a width attached to an imagined absorber surface at a height I. Here a denotes the 
rms height roughness of the absorber. Therefore, we give the loss rate of the n th bound state in terms 
of the general description of scattering processes as a function of the probability of neutrons to dwell 
within the roughness surface region of the absorber/scatterer as 

r„(0 = aioss.n ■ / dz\ip n (z)\ 2 . (9) 

Here we used the fact that the geometry of the wave guide system with and without gravity allows 
for coordinates {x, y, z) with z denoting the transverse coordinate for which the Schrodinger equation 



becomes separable with the product ansatz 



4>n{r) = ip n , g (z) ■ cp xy (x,y) (10) 
with:cf )xv (x,y) = —L=-e tk * x+lk yy . 



Since losses due to non-specular scattering should only depend on local quantities of the surface and 
the probability of finding neutrons at the surface, we assume that the micro-physics of neutron loss 
is independent of the 'macro'-physics of the wave function behavior. The above ansatz represents this 
fact in its product structure which separates the micro-physical quantity ai OSSin completely from the 
behavior of the wave function described by the probability density |?/>n(z)| 2 . Therefore, it is the specific 
loss rate an OS s,n which would depend in a micro-physical model on the roughness properties a (roughness 
variance) and £ (roughness correlation length), i.e, ai OS s,n = ccioss,n(c, £)■ Furthermore this argument 
requires that micro-physical quantities should not depend on the neutron state number. Thus, we have 
aioss,n = aiossVn. Thus, we have 

r„(0 = ai oss • f dz\ip n (z)\ 2 . (11) 

Jl-2a 

We fit this quantity ai OS s to the data. 

This leaves us with the task to determine the integrals f, 2 dz\ip n (z)\ 2 for a state n in a given 
experimental setup. 



A. The wave guide system with gravity 



The bound states i/) n ,g in the linear gravitational potential are confined by two very high potential 
steps above (absorber) and below (mirror). They can be given analytically as: 

Vv 3 M = A n ■ Ai (z/R - C„(0) + B n ■ Bi (z/R - £„(/)) (12) 

and they fulfill the boundary conditions: 



=0 



= A 



= (13) 



which account for the high potential steps bounding the potential from below and above. Eq.s i|12|) and 
(I13fl together determine the energy eigenvalues Cn(0 as functions of the absorber height / which is done 
numerically. 

For a given state ip n ,g the Ai-part in the wave function is exponentially decaying for z > z n (l) = 
£ n (l)R while the Si-part grows exponentially in this region. At z = I both parts add up to zero. As 
long as I > z n (l) we have B n « A n because the Si-part has to compensate only for the exponentially 
small value of the Ai-part at z = I. This, however, implies that at z < I we have \B n Bi(z/R— £n(0)l < - 
A n Ai(z J R — Cn(l))- Therefore, we can ignore the £?i-part of the wave-function in the calculation of an 
asymptotic expression for the r„ jS (Z) at large I. As a result, we have 

^ n , g (z)^A n -Ai(z/R-Cn{l)) (14) 

at / 3> z n . The main effect of the i?i-part is the increase of Cn(0 > z n aX I < z n where the z n are given 
by Eq.El Here we denote all quantities at I ~ » oo with £„, A n etc. while the corresponding quantities of 
the realistic states Eq. Elat finite I are denoted with ( n (l),A n (l) etc. 

The states VVs are we ^ approximated in the large I regime by the asymptotic WKB states 



^, 9 -r4=-(C-C„)- 1/4 -e-i«-« 3/2 
Z\/ttR 



if C>Cn 



(15) 



In the limit I — > oo we have the C«(0 ~~ ¥ Cn given by Eq. EJto 1 st order in WKB and A n = irBi(—( n ). 
(For I — > oo we have ipn,g{z) = A n ■ Ai(z/ R — £»»)• Then the quoted value of A n which is exact ensures 
that (il> n ,g\il>n,g) = 1-) In the case I = R(i < z n consider the WKB-expression for the energy eigenvalues 



dCVCn(0-C = 7r(n-l/4) . (16) 



With 

<l<Cn(l) 



(17) 



we arrive at 



7r 2 (n-l/4) 2 



Cn(0 = r 2 (18) 

which coincides with the known box state expression for large n. 

Now plugging the asymptotics of the gravitationally bound states Ea.lTSlmto Eq.^Jone arrives at 
a prediction for the loss rate that reads 



r„, 9 (7) = aioss • / dz ■ \ip n>g (z)\ 



2 



l-2a 

«i°ss ■ ' \ A n? ■ € * R == for Z»z„. (19) 
2tt# ^/(Z - z „)/i? 



Now look at the behavior of the total neutron flux F through the wave guide if a large number of 
states contribute to it (semi-classical limit). We have F = J A jx ~ 2 n (V , n|V'n) where we get (V'nIV'n) 
from Eq. |S] and A denotes the wave guide cross section. If gravity is present the behavior ri oss ~ 
exp(— 4/3 • (C — Cn(0) 3 ^ 2 ) leads to the fact that each time when I ~ z n = R( n a new state rapidly starts 
to contribute to the transmission. Therefore at a given large height I the number of states contributing 
to <p(l) reads from Eq. [fj] 

I = z n ~ N 2/3 AT(i) - Z 3/2 

which yields asymptotically (if N(l) large and thus AiVv_(jv_i)/iV = 1/JV — > 0) the classical behavior 
in a gravitational field. A simple phase space argument |2| shows that a perfect absorber at the top in 
presence of a gravitational field yields a classical transmission 

F(l)~l 3 / 2 . 



B. A wave-guide system without gravity 



Now we can use the model to derive the loss rate for bound states in the absence of gravity. These 
bound states are well approximated by those which describe the quantum dynamics of a particle in a 
one-dimensional box with infinitely high walls, i.e. the so-called box states. They are given by 

/ 2 . mr 

4>n = yysm(— -z) (20) 

which yields a loss rate given by 

T n (l) = a loss - / dz ■ \-4>„(z)\ 2 

Jl-2a 

'2a 1 . f Ann \\ 

= "loss • —, — sin — — - a) . (21) 

I 2nn \ I J 

A comparison of the full box state expression (line 2 in the above equation) with the numerical result 
for the loss rate with gravity using the full states Eq. ^|is given graphically in Fig. EI 

The exact expression in the second line of Eq. 1211 approaches a constant for n — > oo rendering 
the sum in Ea. 1231 divergent. Therefore, in a realistic fit we have to include the fact that in any real 
experiment the number of box states N in the wave guide is finite. Firstly the collimator system in 
front of the wave guide yields an input vertical velocity distribution of finite width. Secondly for wider 
vertical velocity spectra all neutrons with v z > v^ rlt — 4.3 • m ■ s _1 , the critical velocity of glass, will 
enter the mirror or the absorber directly without forming bound states in the wave guide. In both cases 
the number of box states which is populated by the entering flux of neutrons behaves like N ~ I. If the 
collimator is tuned to yield an input vertical velocity distribution of small width (i.e. about 10 cm/ s) we 



have N/l ~ 2 • [imT 1 corresponding to about N 200 box states in a wave guide of I ~ 100 /im width. 
Thus, we have evaluated a finite sum with N/l — 2 • [imT 1 when comparing the gravity- free prediction of 
Eq. [231 with the experimental data sets (for I < 100 /im we have N < 200 box states populated with N 
approaching w 200 for I — > 100 /im). The dependence of our result on the choice of the cutoff N is very 
weak. If the critical velocity of glass defines the cutoff this results in N/l ~ 20 • /im" 1 (corresponding to 
about N m 2 ■ 10 3 box states in a wave guide of I ~ 100 /im width). We find that the predictions agree 
for both cutoff choices with each other to far better than the experimental accuracy within the Grange 
of the measurement. 

Now in the absence of gravity the asymptotic behavior of the transmitted flux F(l) carried by the 
box states at large I can be given directly from the Ea.s Inland 1211 For na/l <C 1 we can approximate 
Eq. El with 

r„(Z) = ai oss • 16^ 2 /3 • a 3 • ^ if: ^ <C 1 . (22) 
If the approximation in Eg. 1221 were valid for all n we would have 
F(l) ~ ^ e-rnCO-ifugM 

n>l 

= e"^-™ 2 , with : 7 = 16tt 2 /3 • a 3 ■ l~ 3 ■ t mght 

n>l 

= 1.(03(0,0-1) 

- l 3/2 for large/ . (23) 

Here i?„ (g, u) denotes the elliptic theta function where we used Mathematica to evaluate the sum. 
This result differs from the naive classical behavior of a gravity- free wave guide with a perfect absorber: 
In the case of the linear trajectories describing classical particles in absence of the gravitational field, 
we find [H 

F(l) ~ I 2 (24) 

which is easy to imagine since one factor of I obviously has its origin in the relation F(l) ~ A ~ I while 
the 2 nd factor encodes that the range of vertical velocities Av z of particles which pass the wave guide 
without ever touching the absorber also behaves like Av z ~ I. 

Finally, from this situation we expect in general an interpolating behavior of the gravity-free trans- 
mission rate with respect to its power-law dependence on the absorber height. In fact, if the gravity- 
free prediction of Eq. [33] is carried out using the exact expression for the loss rates T n in Eq. for 
N/l > 200 • /im -1 , i.e. N > 2 • 10 4 box states at I w 100 /im, we find that the transmission begins to 
deviate from an Z 3 / 2 -power law towards an /"-dependence with n — > 2 which is the general dependence 
to be expected both classically and quantum mechanically. Thus, the behavior of the gravity-free pre- 
diction as F(l) ~ l 3 t 2 in our given experimental situation is an artifact caused by the relatively small 
number of box states (N < 2 • 10 3 for I < 100 • /im) in the wave guide. 



C. Reversed geometry 



We now turn to the third case of g — ► — g instead of g — * 0. This inversion of gravity is equivalent to 
a setup geometry where the absorber/scatterer is placed at the bottom at z = and a movable mirror 
at z = I above the absorber. For this situation we can follow the derivation of Subsection IIV Al Since 
the absorber is now at z = we have to evaluate the probability integral at this position which implies 
for large / the use of the asymptotic WKB expression 



This results in 



(Cn " 



-1/4 



sm 



if C<Cn 



(25) 



i(n,S,«v.) 



(0 = Q'lo 



2ttR 



1 n| 3 V R R 2 



(26) 



Note, that this loss rate of the reversed geometry is practically independent on the state number n 
since lim n _ 1 .( >0 |-A n | 2 -\/Ci = tt and even at n = 1 it is Pa/Ci/tt — 1 0.5%. Comparing this result 
with the corresponding expression Eq. 1191 for the normal geometry we find that the ratio of the fluxes 

F„ Tev '(l)/Fb 9 \l) of the n th state for I > z n between the reversed and the normal geometry is given 

by 



Fk s> (l) 

_ _ _L 

2ttR ' l^ 1 1 3 V R ' W ' v h 



Since ■ \ Ai\ 2 ■ ^-v/lF' W ' ~ ^ ' 10 5 s for experimental values of L = 0.13 -m, v hor , w 10 -m/s 



and a = 0.75 fim a value of ai oss O(10 4 s _1 ) fitted from a measurement of the normal geometry would 
result in a huge asymmetry under a 7r-rotation around the optical axis of the wave guide when comparing 
the normal and the reversed geometry. 



V. A FIT TO THE DATA 



For a comparison with the measurements we plug the resulting loss rates into the general prediction 
for the transmitted neutron flux Eq. |2HI Together with the repopulation coefficient pi accounting for 
eventual shifts between split bottom mirrors one predicts 



F(l) = F Q + J2 Fn(D = F Q + C • <L • e-r&W-^ + £ . e - r >-»-^ I (28 ) 

ra I n>l ) 

Fq is the detector background and C the total flux normalization. ai oss is the universal parameter 
introduced above which parameterizes the scatterer strength. This quantity is in general expected to 
be a rather weak function of the roughness parameters a and £ which would be determined in principle 
by a microphysical calculation of the loss due to non-specular scattering at the rough absorber surface. 
Thus it does not depend on the absence or presence of gravity nor on the state number n. 

In the case of I < z n , the absorber/scatterer will begin to squeeze the bound states once they start 
to 'feel' it sufficiently strongly. This implies further that due to E ntg (l) > E n , g (oo) — E n<g pure grav ity 
for the true bound states a sufficiently small I leads to E n> g(l) mg ■ I. Therefore, the calculation of 

the r["g S is done using the full realistic bound states Ea. 1121 and deriving the relations corresponding to 
Ea.s ll9l 1211 and numerically which incorporates the mentioned behavior. 



A. Normal geometry 



The detector background 0o has been measured independently to yield Fq = (0.0043 ± 0.0004)s _1 
for the 1999 measurement 0] and Fq = (0.0004 ± 0.0001)s _1 for the new 2002 measurement using an 
improved setup |lfij |. Thus, one remains with having to determine the two universal quantities ai OS s and 
C from the data. C turns out to be completely fixed by the data points at I > 70 /im, and thus also 
has been measured. The fit therefore will be a 1-parameter one determining ai oss as long as all the 
populations p n stay to be equal. 

We fit now Eq. EHI to the newer data from the run of the experiment in 2002 which has a different 
absorber and better statistics, and systematic effects are smaller than in 1999. The results are shown 
in Fig. 2. The value of ai OS s is found in a fit to the data. The result of the fit yields (L — 13 cm, 
Vhor. ~ 5 • m - s _1 ) 

ai oss = {3A±0.1)-I0\s- 1 . (29) 

The fit was done using neutrons with only one value of the horizontal velocity which was chosen to 
be the average velocity v^or. ~ 5 • m • s _1 . This approximation produces essentially the same results 
as if one uses the full actual spectrum of horizontal velocities produced by the collimating system to 
calculate the transmission. 

The 2002 run was performed with only one bottom mirror so that no significant repopulation effects 
of the ground state are expected. However, it was necessary to allow the population of the ground state 



to shift towards 75% compared to the excited states in order to describe the data: p\ = 0.77 ± 0.04. 
This suppression is significant, but it is not large leaving all the states still to be approximately equally 
populated. 

In the measurement of the 1999 run, two bottom mirrors were used. It was tried to shift these 
mirrors relative to each other vertically and in alignment by a few micrometers. Small shifts of a few 
microns between two bottom mirrors change the population of the ground state and the next state quite 
drastically: A shifted geometry with a relative mirror shift of about 5 /im and no relative tilt of the two 
bottom mirrors results in p\ 0.25 compared to p n — 1, n > 2. The fit to the data (see Fig. results 
in aioss = (5.3 ± 0.5) • 10 4 s _1 and p\ — 0.24 ± 0.1 which would be consistent with the possible relative 
mirror shift discussed above. 



B. Reversed geometry 



In a second setup described by the term 'reversed geometry', the absorber was placed at the bottom 
and the mirror above. Here, the position where the scattering-inducing roughness is found at z = 0. 
This case can also be derived from the states Eq. 1121 by just placing the roughness appropriately on 
bottom and then calculating the corresponding loss rates. The bound states are now again given by the 
Airy functions. However, at z = where the absorber is now placed, they do not decay exponentially 
fast and the gravitationally bound states will be strongly absorbed at arbitrary heights I of the mirror 
at top - quite contrary to the normal setup. If gravity were absent, such a 7r-turn of the wave guide 
around its optical axis would not have such an effect. Thus, there would be no preferred direction in 
space forcing the transmission factor to be invariant under rotations around the optical axis of the 
system since the box states describing the simple confinement situation without gravity are symmetric 
with respect to the optical axis. This check was done with the 1999 data. We took the prediction for 
gravitationally bound states of Eq. 1121 in the normal setup and fitted ai oss and p\ to the data. The 



model yields p\ = 0.24 ±0.1 and ct\ 



(5.3 ±0.5) 



10 4 £ 



Taking now these values for ai oss and 



pi from the fit, we can calculate a prediction for the setup with reversed geometry using the results of 
Subsection IIV CI - which is now entirely fixed and not fitted any more. Using the experimental setup 
parameters used there, L = 13 cm and «hor. ~ 10 • m/s we can estimate the suppression factor using 
Eq. E3 on the ground state which yields 



V) 



0.03 



(30) 



This fits well with the two data points obtained experimentally in this reversed setup. Using the full 
prediction calculated again numerically from Ea. ll2l in the reversed setup, the comparison with the 1999 
data is shown in Fig. 01 In treating the 1999 data other models 0, El showed more pronounced steps. 
These steps are not present in the (more advanced) model presented in this paper. 
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FIG. 3: Behavior of the loss rate V (plotted in ar- 
bitrary units) for the pure box state n = 5 (dash), 
an approximate power law, and for the gravitationally 
bound state n = 5 found numerically from using the 
full states Ea. 1121 fsolidi. showing exponentially fast 
decay above some I = O(R). 
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FIG. 4: Open Circles: Transmission as a function of the 
absorber-mirror distance. Filled Circles: Transmission 
measured with inverse geometry - the absorber on the 
bottom and a 10 cm mirror on top (both are 1999 
measurement data @|). 



The above described asymmetry of the neutron transmission of a wave guide with one absorbing 



and one reflecting wall under rotations around its optical axis in the presence of gravity has indeed 
been measured and can be described quantitatively. This rules out the possibility to explain the data 
just by confinement effects in terms of box states inside a rectangular box-shaped potential 0] and 
thus establishes the gravitational nature of force which binds the neutron bound states to the bottom 
mirror. For further discussion, see also |2(1 |. The measurement of the reverse geometry is a very good 
test of our (or probably any) model, since the result depends strongly on the absorber properties. The 
fact that the absorber model which describes the data in the normal geometry can also be used for the 
reversed geometry, without any new fit parameters, gives us confidence that our model is a useful tool 
to describe the data. 



VI. SUMMARY 



For the first time the existence of quantized bound states of neutrons in the gravitational field of 
the earth above a horizontal glass mirror was experimentally demonstrated. Here we present a quantum 
mechanical model providing an accurate fit to the data. The difficult part is the incorporation of the 
neutron absorber into this model. We show that we can describe its action with only one phenomeno- 
logical fit parameter in several configurations. The most important configuration has bound quantum 
states in the linear gravitational potential. We consider repopulation effects of these states when a step 
from a second mirror is present. The emerging reduction of the first state has been calculated. The 
standard setup uses only one mirror. Also in this case, the data show an unexpected reduction of the 
first quantum state as well. 

Further, it is shown that the experiment would generate different results, if gravity were hypothet- 
ically to be turned off. This configuration is a wave-guide system with box states which describe the 
quantized motion of a particle a one-dimensional box with infinitely high walls. 

A striking difference is found, if the gravitational field changes sign. This last case has been observed 
by our measurement and can be described quantitatively within the same modeling of the absorber 
mechanism presented here. Thus, the dominating effect of gravity on the formation of the bound states 
has been demonstrated since the measurement did not show just a simple confinement effect, as proposed 
in 0| . On the contrary, it could not be described theoretically, if the earth's gravitational field was not 
present. 
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